Large Bychkov-Rashba spin-orbit coupling in high mobility GaN/AlGaN 

heterostructures 
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We present low temperature magneto-conductivity measurements of a density-tunable and high 
mobility two-dimensional electron gas confined in the wide bandgap GaN/AlGaN system. We ob- 
served pronounced anti-localization minima in the low-field conductivity, indicating the presence 
of strong spin-orbit coupling. Density dependent measurements of magnetoconductivity indicate 
that the coupling is mainly due to the Bychkov-Rashba mechanism. In addition, we have derived 
a closed-form expression for the magnetoconductivity, allowing us to extract reliable transport pa- 
rameters for our devices. The Rashba spin-orbit coupling constant is a so ~ 6x 10" 13 eVm, while 
the conduction band spin-orbit splitting energy amounts to A so ~ 0.3meV at n e = lx 10 16 m -2 . 



GaN has emerged as a leading material for a vari- 
ety of new device applications, ranging from solid-state, 
ultra-violet optical sources to high power electronics Q. 
While the performance of many devices fabricated from 
GaN has been stunning, several fundamental physical 
processes remain to be understood. A prime example is 
spin-orbit coupling in GaN and its heterostructures. The 
burgeoning field of spintronics has invigorated the study 
of spin-orbit coupling in semiconducting materials |2j, yj . 
To date, much experimental effort has been devoted to 
narrow bandgap material systems (InAs, InGaAs, GaAs, 
etc) as spin-orbit coupling is expected to be strong in 
these systems. Conversely, far less experimental effort 
has been directed toward wide bandgap systems like GaN 
in which spin-orbit effects are predicted to be suppressed 
by the large fundamental bandgap, E g , and reduced spin- 
orbit splitting, Ao, of the valence band at zone center. 
Indeed, in the kp formalism the bare Rashba spin- 
orbit coupling constant for electrons, ao, scales as: «o~ 
Ao/E 2 ,. As the value of Ao for GaAs exceeds that of 
GaN by a factor of 30, it is reasonable to suspect that 
spin splitting of the conduction band in GaN-based het- 
erostructures would be insignificant compared to GaAs 
and other narrow gap heterostructures. 

Spin-orbit coupling for conduction band electrons in 
bulk GaN was considered by Krishnamurthy [5| who cal- 
culated that the spin relaxation times in bulk GaN should 
exceed the spin relaxation time in GaAs by three orders 
of magnitude, thus making GaN an excellent candidate 
for transport of spin polarized currents over macroscopic 
distances. However, in Ref. []| GaN was assumed to 
have the zinc-blende lattice structure. GaN is typically 
grown in the more stable wurtzite phase. It is known that 
the symmetry of the underlying crystal has a profound 
impact on spin-orbit induced splittings in the conduc- 
tion band 0,0 ■ While the work of Ref. [j| is suggestive, 
very few experimental results for bulk wurtzite GaN have 
been reported || and the impact of spin-orbit coupling 
on transport in wurtzite GaN/AlGaN heterostructures 
remains an open question. A few preliminary experi- 



ments have considered spin-orbit cou pling for the two- 
dimensional electron gas (2DEG) 0, lid fill fl2T | in a 
narrow parameter space of high density and low mobil- 
ity. The influence of spin-orbit coupling in the limit of 
low 2DEG density and high mobility has not been ad- 
dressed. Furthermore, probative experiments in which 
conductivity is tuned over a broad range, and theoretical 
analysis specifically tailored to the physics of GaN are 
needed to understand the mechanisms of spin-orbit cou- 
pling and accurately quantify spin-orbit effects in high 
mobility GaN 2DEGs. 

In this Letter we present an analysis of low tem- 
perature magnetoconductivity measurements in a series 
of high mobility 2DEGs confined in the wide bandgap 
GaN/AlGaN system. Experiments are conducted with 
gated Hall bars that allow access to a previously in- 
accessable range of low density, 5 x 10 15 m~ 2 < n e < 
1.8 x 10 16 m~ 2 , and very high mobilities 1.4m 2 /Vs < 
fi < 8.7m 2 /Vs. We observe non-monotonic behavior in 
the magnetoconductivity with a pronounced antilocal- 
ization minimum at B~ 2mT, indicating the presence 
of significant spin-orbit coupling. In addition, we have 
derived an exact closed-form expression for the magne- 
toconductivity, that allows for the extraction of reliable 
spin-orbit parameters relevant to our devices. The rela- 
tive simplicity of the formula for the magnetoconductiv- 
ity greatly facilitates data fitting. Importantly, the mag- 
netic field at which the magnetoconductivity minimum 
occurs does not depend sensitively on electron density. 
As we shall show, this result implies that the Bychkov- 
Rashba mechanism is the dominant spin-orbit coupling 
in our samples. The extracted Rashba coupling constant 
a so =6x 10 _13 eVm is large, resulting in spin-split ener- 
gies ranging from 0.2meV to 0.4meV within the density 
range of our experiment. The value of a so in GaN is 
comparable to that seen in the narrower bandgap GaAs 
system. Our findings place severe constraints on the use 
of GaN heterostructures for polarized spin transport, but 
also suggest that GaN may be implemented in applica- 
tions where only narrow bandgap materials have been 
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considered previously. 

Magnetoconductivity in two dimensions has been stud- 
ied extensively in the diffusive limit [T3L fl4| . The con- 
ductivity of a 2DEG in classically weak magnetic fields, 
ct(-B), shows signatures of quantum interference that de- 
pend on the magnetic field and spin-orbit coupling. Spin 
relaxation due to spin-orbit coupling and impurity scat- 
tering produces a positive contribution to the conduc- 
tivity known as antilocalization. Magnetic field sup- 
presses this antilocalization. The functional dependence 
of Aa(B), where Aa(B) = <r(B) - a{B = 0), at small 
magnetic fields depends on the relative contributions of 
dephasing (characterized by a dephasing rate 1/t v ) and 
the spin-relaxation rate given by 1/ t so . Strong spin relax- 
ation compared to dephasing leads to a pronounced min- 
imum in the magnetoconductivity for excursions away 
from B=0. The observation of an antilocalization mini- 
mum in the magnetoconductivity is the signature of spin- 
orbit coupling. 

The Bychkov-Rashba interaction in a 2 D sy stem can 
be described by the following Hamiltonian 



j sy s 



H = h a so a ■ (z x p). 

Am 



(1) 



Here, a so is the spin-orbit coupling strength, the split- 



ting energy is A so =2a so pF (at the Fermi surface). In 
a 2D system with spin-orbit interactions, the dominant 
spin relaxation process is typically the D'yakonov-Perel' 
(DP) mechanism ^||. It describes the relaxation of the 
electron spin in the presence of a spin-splitting, A so . Re- 
laxation occurs because the direction of the axis of spin 
precession is tied to the direction of the electron momen- 
tum, which changes randomly with each collision. As a 
result, the net precession after N collisions is typically 
\fN Asot/% in the diffusive transport regime, where t 
is the mean free time. Consequently, the time it takes 
to randomize the spin is t so ~ h 2 / A 2 so t. Expressed in 
terms of a so , l/T so =D(2ma so /H) , where D = v f t/2 is 
the diffusion constant in 2D. 



We analyze the measured magnetoconductivity data 
with an analytical formula we have derived for Aa(B) in 
the presence of the Bychkov-Rashba interaction that is 
valid in the diffusive regime. The details of our derivation 
are presented in an upcoming publication |l7j . This for- 
mulation greatly simplifies extraction of the transport pa- 
rameters from experimental data. The relevant magnetic 
field scales are: B so = H/4eDT so , B v = fi/4e£>T ¥ ,, and 
Btr=ti/4£DT. The diffusive limit is defined as B -C B tr . 
In this limit, Aa(B) is independent of B tr and the ex- 
pression for A<t(B) reads: 



Aa(B) 



C 
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The values of u s and v K are: 
= 2 5 cos 



27T, 
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where the variables 5 and 9 are equal to: 



1 - Absob? - b 2 so 
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b so 



(4) 
(5) 
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ip(z) is the di-gamma function, while b so =B so / B and 
b v =B v /B. The constant C is such that Aer(0) = 
0. Eq. (2) provides an analytical, closed-form, solu- 
tion for the magnetoconductivity in the presence of the 



Bychkov-Rashba interaction. The result is expressed in 
terms of the two parameters B so and B v . Because, 
B so = (ma so ) 2 j eh, determining B so as a fitting parame- 
ter gives, if the mass is known, directly the value of the 
spin-orbit coupling a so defined in Eq. (1). 

The samples used in this study arc single interface 
GaN/AlGaN heterostructures grown by plasma- assisted 
molecular beam epitaxy (MBE) on (0001) oriented GaN 
templates. After an initial 1/mi GaN buffer layer, a 16nm 
thick Al 2 ,Gai_3;N barrier layer (x varies between 0.08 and 
0.12) is grown followed by a 3nm thick GaN capping 
layer. The 2DEG is formed at the lower GaN/AlGaN 
interface without modulation doping due to the effects 
of spontaneous and piezoelectric polarization 0|- Hall 
bars with 100/mi width and 2mm length are defined with 
14 voltage probes symmetrically placed along the device. 
A Ni/Au gate, separated by 50nm of Si02 from the wafer 
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FIG. 1: Transport parameters of our gated Hall bar devices. 
Sample A (solid black circles) has lower conductivity, placing 
it in the diffusive transport regime. Sample B (solid black 
squares) is a high mobility device. The combination of the 
two devices allows access to a broad range of 2D conductivity. 

surface, is used to control the 2DEG density. The mag- 
netic field dependent conductivity was measured with 
standard low frequency lock-in techniques at T=0.3K. All 
investigated samples show an unambiguous minimum in 
low field magnetoconductivity, indicative of strong spin- 
orbit coupling. 

We discuss in detail two samples with different con- 
ductivity (see Fig. 0}. The first structure, sample A, is 
designed to have a mobility and carrier density which 
place it in the diffusive limit at low carrier density. With 
Sample A, we are able to tune the conductivity from the 
diffusive limit, where our theory is strictly valid and the 
extracted parameters are most accurate, to the ballistic 
regime at higher carrier density. Our objective is to ex- 
tract reliable parameters for sample A in the diffusive 
limit and then monitor the evolution of the magnetocon- 
ductivity as the transport moves into the ballistic regime. 
Sample B, which has fi = 8.7m 2 /Vs at n e =1.8 x 10 16 m~ 2 , 
is a very high mobility sample [2lL |22j which extends our 
study deep into the ballistic regime. 

The experimentally obtained data are plotted as a 
change in conductivity Aa(B) = a(B) — a(B = 0) and 
are shown in Fig. [2] for sample A. We note that Eq. (2) 
predicts a crossover field from negative magnetoconduc- 
tivity to positive magnetoconductivity at B f=a B so . Sev- 
eral features are evident in the data. While we change the 
density by a factor of two and the conductivity by a factor 
of four, the field scale at which the magnetoconductivity 
minimum occurs does not change within our experimen- 
tal accuracy. This fact has two immediate consequences. 
As the magnitude of the spin-orbit coupling constant a so 
is directly proportional to B SOl the data implies that the 
Bychkov-Rashba coupling strength does not change as a 
function of density. In addition, the lack of field depen- 
dence of the conductivity minimum justifies our neglect 
of the linear and cubic Dresselhaus terms in the model 
Hamiltonian, since the presence of Dresselhaus coupling 
dictates that the field scale B so acquires an explicit den- 
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FIG. 2: Measured magnetoconductivity data (symbols) of 
sample A for small magnetic fields over a wide range of den- 
sity and mobility (inset). In all cases Aa(B = 0) = 0. The 
traces have been shifted for clearer presentation. The solid 
lines denote the fit of Eq. (2) at each density. The verti- 
cal bars are at the position of B tr , denoting the limit to the 
diffusive regime. 



sity dependence, B so ~ 7 2 n 2 [XU EI3| • 

We also see that while the magnetic field position of the 
conductivity minimum does not change with conductiv- 
ity, the amplitude of the antilocalization effect docs. The 
antilocalization induced drop in conductivity is largest at 
highest conductivity while it is suppressed at lower con- 
ductivity. Since the field scale of the conductivity mini- 
mum, B so , is largely constant, the reduction in the am- 
plitude of the antilocalization behavior can be attributed 
to the change in B v as the density and mobility are re- 
duced. 

Using the model presented earlier, the changes to the 
conductivity due to quantum interference in the pres- 
ence of spin-orbit coupling at each density were calcu- 
lated (solid lines in Fig.[2J), allowing for the extraction of 
B so and B v . Also shown by vertical lines are the values 
of B tr for each density. As expected, the deviation of the 
fits from the data becomes significant at fields B ~ B tr . 

The values for B so and B v extracted from the magne- 
toconductivity are plotted as a function of 2DEG density 
in Fig. The field scale B v associated with dephasing 
decreases its value as the conductivity increases. The 
values for B v correspond to a dephasing time r^lOOps 
in sample A. As could be expected from inspection of 
the raw data, B so remains at ~2mT for the whole range 
of densities explored. The corresponding values for the 
spin dephasing time t so vary from 2ps (high density) to 
lOps (low density). Using B so ^2mT, the spin-orbit cou- 
pling parameter cv so =6x lCP 13 eVm is calculated, yielding 
a density dependent splitting, A so , of the two spin sub- 
bands at the Fermi edge of 0.2meV - 0.3meV for this 
sample. 
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FIG. 3: The density dependent values of B so and B v ex- 
tracted from the fits in Fig. [21 While B v decreases by a factor 
of six with increasing density, B ao shows little density depen- 
dence. The inset shows the evolution of r so as a function of 
2D density. 
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structures is nearly identical. Given a maximum density 
n e =1.8x 10 16 m -2 , a zero field spin splitting of 0.4meV is 
calculated for sample B. It is interesting to compare our 
value of spin splitting to the results of Chou et al. j2^| 
who examined the zero-field splitting observed in quan- 
tum point contacts (QPC) fabricated on a similar high 
mobility GaN heterostructure. In a QPC built on a GaN 
2DEG with n e = lxl0 16 ur 2 and /i=5.6m 2 /Vs, Chou ob- 
served a zero-field splitting of 0.39meV, a result which 
compares favorably with our value determined from mag- 
netoconductivity measurements of a GaN 2DEG. 

The origin of the large Bychkov-Rashba coupling in 
GaN heterostructures is yet to be fully explained. One 
possible explanation is found in the extremely large elec- 
tric fields at the AlGaN/GaN interface generated by the 
polarization discontinuity. Self-consistent calculations 
suggest that the electric field in our GaN heterostruc- 
ture with n e =5x 10 15 m -2 , is approximately 10 times the 
value in an equivalent density GaAs 2DEG structure [24| • 
Since a so =aoE, where ao is the fundamental spin-orbit 
coupling parameter for a particular material system and 
E is the value of the electric field at the heterointerfacc, 
it may be that large zero-field spin splitting is due to 
the large increase in electric field. It is also possible that 
the fundamental spin-orbit coupling may be enhanced in 
wurtzite GaN due to coupling between the conduction 
band and higher energy conduction bands. A calcula- 
tion of ao in GaN which accounts for the wurtzite sym- 
metry and includes remote band effects is still lacking. 
While further study is needed, our experiments clearly 
indicate that spin-orbit coupling dramatically influences 
the transport properties of high mobility GaN 2DEGs. 

* Corresponding author: manfra@lucent.com 



FIG. 4: Magnetoconductivity data for sample B. The high 
mobility of sample B dictates that the transport is ballistic. 
As with sample A, the conductivity minimum is situated near 
B=2mT and does not depend significantly on 2D density. 



We now turn to sample B. In sample B, B tr covers 
the range between lOO^jT and 440/xT and is thus signifi- 
cantly smaller than the field scale B so . It is clear from the 
forgoing discussion that our closed form expression for 
the magnetoconductivity has limited applicability in this 
regime. We therefore do not attempt to fit the data. Nev- 
ertheless, we still can glean important information from 
the raw data. The magnetoconductivity data for sample 
B are presented in Fig. The conductivity minimum 
does not change position with density and the minimum 
conductivity is again near B=2mT, as in sample A. It 
follows that the value of a so for sample B will be approx- 
imately the same as in sample A. Indeed, it is not supris- 
ing that B so has not changed significantly in moving from 
sample A to sample B as the layer sequence of the two 
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